Plastid phylogenomics and fossil evidence provide new insights into the evolutionary complexity of the 'woody clade' in Saxifragales.

BACKGROUND: The "woody clade" in Saxifragales (WCS), encompassing four woody families (Altingiaceae, Cercidiphyllaceae, Daphniphyllaceae, and Hamamelidaceae), is a phylogenetically recalcitrant node in the angiosperm tree of life, as the interfamilial relationships of the WCS remain contentious. Based on a comprehensive sampling of WCS genera, this study aims to recover a robust maternal backbone phylogeny of the WCS by analyzing plastid genome (plastome) sequence data using Bayesian inference (BI), maximum likelihood (ML), and maximum parsimony (MP) methods, and to explore the possible causes of the phylogenetic recalcitrance with respect to deep relationships within the WCS, in combination with molecular and fossil evidence. RESULTS: Although the four WCS families were identically resolved as monophyletic, the MP analysis recovered different tree topologies for the relationships among Altingiaceae, Cercidiphyllaceae, and Daphniphyllaceae from the ML and BI phylogenies. The fossil-calibrated plastome phylogeny showed that the WCS underwent a rapid divergence of crown groups in the early Cretaceous (between 104.79 and 100.23 Ma), leading to the origin of the stem lineage ancestors of Altingiaceae, Cercidiphyllaceae, Daphniphyllaceae, and Hamamelidaceae within a very short time span (∼4.56 Ma). Compared with the tree topology recovered in a previous study based on nuclear genome data, cytonuclear discordance regarding the interfamilial relationships of the WCS was detected. CONCLUSIONS: Molecular and fossil evidence imply that the early divergence of the WCS might have experienced radiative diversification of crown groups, extensive extinctions at the genus and species levels around the Cretaceous/Paleocene boundary, and ancient hybridization. Such evolutionarily complex events may introduce biases in topological estimations within the WCS due to incomplete lineage sorting, cytonuclear discordance, and long-branch attraction, potentially impacting the accurate reconstruction of deep relationships.

The Complete Mitochondrial Genome of Paeonia lactiflora Pall. (Saxifragales: Paeoniaceae): Evidence of Gene Transfer from Chloroplast to Mitochondrial Genome.

Paeonia lactiflora (P. lactiflora), a perennial plant renowned for its medicinal roots, provides a unique case for studying the phylogenetic relationships of species based on organelle genomes, as well as the transference of DNA across organelle genomes. In order to investigate this matter, we sequenced and characterized the mitochondrial genome (mitogenome) of P. lactiflora. Similar to the chloroplast genome (cpgenome), the mitogenome of P. lactiflora extends across 181,688 base pairs (bp). Its unique quadripartite structure results from a pair of extensive inverted repeats, each measuring 25,680 bp in length. The annotated mitogenome includes 27 protein-coding genes, 37 tRNAs, 8 rRNAs, and two pseudogenes (rpl5, rpl16). Phylogenetic analysis was performed to identify phylogenetic trees consistent with Paeonia species phylogeny in the APG Ⅳ system. Moreover, a total of 12 MTPT events were identified and 32 RNA editing sites were detected during mitogenome analysis of P. lactiflora. Our research successfully compiled and annotated the mitogenome of P. lactiflora. The study provides valuable insights regarding the taxonomic classification and molecular evolution within the Paeoniaceae family.

Codon Usage Analyses Reveal the Evolutionary Patterns among Plastid Genes of Saxifragales at a Larger-Sampling Scale.

Saxifragales is a 15-family order of early-divergent Eudicots with a rich morphological diversity and an ancient rapid radiation. Codon usage bias (CUB) analyses have emerged as an essential tool for understanding the evolutionary dynamics in genes. Thus far, the codon utilization patterns had only been reported in four separate genera within Saxifragales. This study provides a comprehensive assessment of the codon manipulation based on 50 plastid genes, covering 11 constituent families at a larger sampling scale. Our results first showed a high preference for AT bases and AT-ending codons. We then used effective number of codons (ENC) to assess a range of codon bias levels in the plastid genes. We also detected high-informative intrafamilial differences of ENC in three families. Subsequently, parity rule 2 (PR2) plot analyses revealed both family-unique and order-shared bias patterns. Most importantly, the ENC plots and neutrality analyses collectively supported the dominant roles of selection in the CUB of Saxifragales plastid genes. Notably, the phylogenetic affinities inferred by both ML and BI methods were consistent with each other, and they all comprised two primary clades and four subclades. These findings significantly enhance our understanding of the evolutionary processes of the Saxifrage order, and could potentially inspire more CUB analyses at higher taxonomic levels.

Comparative and phylogenetic analyses of the chloroplast genomes of species of Paeoniaceae.

Plants belonging to family Paeoniaceae are not only economically important ornamental plants but also medicinal plants used as an important source of traditional Chinese medicine. Owing to the complex network evolution and polyploidy evolution of this family, its systematics and taxonomy are controversial and require a detailed investigation. In this study, three complete chloroplast genomes of sect. Paeonia, one of the sections of Paeonia, were sequenced and then analysed together with 16 other published chloroplast genomes of Paeoniaceae species. The total lengths of the chloroplast genomes of these species were 152,153-154,405 bp. A total of 82-87 protein-coding genes, 31-40 tRNA genes and 8 rRNA genes were annotated. Bioinformatics analysis revealed 61-74 simple sequence repeats (SSRs) in the chloroplast genomes, most of which have A/T base preference. Codon usage analysis showed that A/U-ending codons were more positive than C/G-ending codons, and a slight bias in codon usage was observed in these species. A comparative analysis of these 19 species of Paeoniaceae was then conducted. Fourteen highly variable regions were selected for species relationship study. Phylogenetic analysis revealed that the species of sect. Paeonia gathered in one branch and then divided into different small branches. P. lactiflora, P. anomala, P. anomala subsp. veitchii and P. mairei clustered together. P. intermedia was related to P. obovata and P. obovata subsp. willmottiae. P. emodi was the sister to all other species in the sect. Paeonia.

Rates of niche and phenotype evolution lag behind diversification in a temperate radiation.

Environmental change can create opportunities for increased rates of lineage diversification, but continued species accumulation has been hypothesized to lead to slowdowns via competitive exclusion and niche partitioning. Such density-dependent models imply tight linkages between diversification and trait evolution, but there are plausible alternative models. Little is known about the association between diversification and key ecological and phenotypic traits at broad phylogenetic and spatial scales. Do trait evolutionary rates coincide with rates of diversification, are there lags among these rates, or is diversification niche-neutral? To address these questions, we combine a deeply sampled phylogeny for a major flowering plant clade-Saxifragales-with phenotype and niche data to examine temporal patterns of evolutionary rates. The considerable phenotypic and habitat diversity of Saxifragales is greatest in temperate biomes. Global expansion of these habitats since the mid-Miocene provided ecological opportunities that, with density-dependent adaptive radiation, should result in simultaneous rate increases for diversification, niche, and phenotype, followed by decreases with habitat saturation. Instead, we find that these rates have significantly different timings, with increases in diversification occurring at the mid-Miocene Climatic Optimum (∼15 Mya), followed by increases in niche and phenotypic evolutionary rates by ∼5 Mya; all rates increase exponentially to the present. We attribute this surprising lack of temporal coincidence to initial niche-neutral diversification followed by ecological and phenotypic divergence coincident with more extreme cold and dry habitats that proliferated into the Pleistocene. A lack of density-dependence contrasts with investigations of other cosmopolitan lineages, suggesting alternative patterns may be common in the diversification of temperate lineages.

Transcriptomic Sequencing and Co-Expression Network Analysis on Key Genes and Pathways Regulating Nitrogen Use Efficiency in Myriophyllum aquaticum.

Massively input and accumulated ammonium is one of the main causes of eutrophication in aquatic ecosystems, which severely deteriorates water quality. Previous studies showed that one of the commonly used macrophytes, Myriophyllum aquaticum, was capable of not only withstanding ammonium of high concentration, but also efficiently assimilating extracellular ammonium to constitutive amino acids and proteins. However, the genetic mechanism regulating such efficient nitrogen metabolism in M. aquaticum is still poorly understood. Therefore, RNA-based analysis was performed in this study to understand the ammonium regulatory mechanism in M. aquaticum in response to various concentrations of ammonium. A total of 7721 genes were differentially expressed, of which those related to nitrogen-transport, assimilation, and remobilization were highly-regulated in response to various concentrations of ammonium. We have also identified transcription factors and protein kinases that were rapidly induced in response to ammonium, which suggests their involvement in ammonium-mediated signalling. Meanwhile, secondary metabolism including phenolics and anthocyanins biosynthesis was also activated in response to various concentrations of ammonium, especially at high ammonium concentrations. These results proposed a complex physiological and genetic regulation network related to nitrogen, carbohydrate, transcription factors, and secondary metabolism for nitrogen use efficiency in M. aquaticum.

Identification and development of microsatellite (SSRs) makers of Exbucklandia (HAMAMELIDACEAE) by high-throughput sequencing.

Hamamelidaceae (Saxifragales, previously Rosales) comprises approximately six subfamily, 30 genera and 140 species, most of which are Tertiary relicts. Exbucklandia is the only genus of the subfamily Exbucklandioideae, Hamelidaceae, containing only 2-4 species. Of them, the species E. longipetala H. T. Chang is endemic to China and listed as endangered in The Biodiversity Red List of China: Higher Plant, yet some taxonomists put forward that E. longipetala should be merged into E. tonkinensis (Lecomte) H. T. Chang. Currently, there was nearly no phylogeographic studies on this genus possibly due to the deficiency of efficient molecular markers. In this study, we sequenced the genome of E. tonkinensis based on high throughput sequencing technology, and obtained approximately 6 G raw data, which was further de novo assembled into 303,481 contigs. Based on them, 15,326 SSRs were identified from 13,596 contigs, and primers were successfully designed for 10,660 SSRs. A total of 139 paired primers were synthesized, 106 of them were successfully amplified in six Exbucklandia individuals with expected PCR product size, and 24 demonstrated to be polymorphic among three Exbucklandia populations. Accordingly, the expected and observed heterozygosity were between 0.097-0.717 and 0.098-0.583. Based on these efforts, future researches on genetic diversity and population structure of Exbucklandia can be performed to understand its phylogenetic origins and phylogeographic pattern.

Resolving the systematic positions of enigmatic taxa: Manipulating the chloroplast genome data of Saxifragales.

Accurately resolving the phylogeny of enigmatic taxa is always a challenge in phylogenetic inference. Such uncertainties could be due to systematic errors or model violations. Here, we provide an example demonstrating how these factors affect the positioning of Paeoniaceae within Saxifragales based on chloroplast genome data. We newly assembled 14 chloroplast genomes from Saxifragales, and by combining these genomes with those of 63 other angiosperms, three datasets were assembled to test different hypotheses proposed by recent studies. These datasets were subjected to maximum parsimony, maximum likelihood and Bayesian analyses with site-homogeneous/heterogeneous models, different data partitioning strategies, and the inclusion/exclusion of weak phylogenetic signals. Three datasets exhibited remarkable heterogeneity among sites and among taxa of Saxifragales. Phylogenetic analyses under homogeneous models or maximum parsimony showed a closer relationship of Paeoniaceae with herbaceous families in the order. Data partitioning strategies did not change the general tree topology. However, PhyloBayes analysis under the CAT+GTR model resulted in a relationship closer to woody families. We conclude that although genomic data significantly increase the phylogenetic resolution of enigmatic taxa with high support, the phylogenetic results inferred from such data might be analysis or signal dependent. The analytical pipeline outlined here combines phylogenomic inference methods with evaluation of lineage-specific rates of substitution, model selection, and assessment of systematic error. These methods would be applicable to resolve similar difficult questions in the tree of life.